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Abstract 

In a dipolar Fermi gas, the dipole-dipole interaction between fermions can be turned into a dipo- 
lar Ising interaction between pseudospins in the presence of an AC electric field. When trapped in 
a 2D optical lattice, such a dipolar Fermi gas has a very rich phase diagram at zero temperature, 
due to the competition between antiferromagnetism and superfluidity. At half filling, the antifer- 
romagnetic state is the favored ground state. The superfluid state appears as the ground state at 
a smaller filling factor. In between there is a phase-separated region. The order parameter of the 
superfluid state can display different symmetries depending on the filling factor and interaction 
strength, including d-wave (d), extended s-wave (xs), or their linear combination (xs + ix d). The 
implication for the current experiment is discussed. 
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I. INTRODUCTION 

The successful creation of 40 K 87 Rb polar molecules with high space density []J has sparked 
a lot of research interests because of its potential applications in various areas |2|. An 
ultracold gas of polar molecules can provide a platform to study strongly-correlated systems 
jsj], due to the confinement flexibility and the control over effective interactions 4|. The 
dipole-dipole interaction can be turned into different effective interactions in the presence 
of a DC or AC electric field. When the dipole direction is fixed by a DC electric field, a 
p-wave superfluid state with the dominant p z symmetry may appear in a 3D dipolar Fermi 
gas [5]. When an AC electric field is applied to couple different rotational levels of molecules, 
a p x + ip y superfluid state 6] and a ferro-electric state [7| may show up in quasi-2D dipolar 
Fermi gases. Creating the antiferromagnetic state and exotic superfluid states such as p- 
wave and d-wave states has been a focused goal in the research with optical lattices. In this 
work we predict that some of these phases can be created in polar molecules trapped in a 
2D square optical lattice in the presence of a resonant AC electric field. 

For the 40 K 87 Rb molecule, rotational excitation energies are much smaller than electronic 
and vibrational excitation energies. In the presence of the resonant AC electric field for the 
first two rotational states, electronic and vibrational excitations can be ignored, and the 
rotational state of the dipole along or against the electric-field direction can represented by 
pseudospin-up or pseudospin-down state. The dipole-dipole interaction is thus effectively a 
dipolar Ising interaction between pseudospins. The strength of this Ising interaction can be 
much larger than the super-exchange interaction in a unpolarized two-component Fermi gas 
trapped in the same optical lattice. In this aspect it is easier to create antiferromagnetic 
and superfluid states with exotic symmetries in the dipolar Fermi gas. 

This work is organized as follows. First we derive the effective Hamiltonian for this 
system. Then the antiferromagnetic state and superfluid states with different symmetries 
are studied in the mean-field approximation. The ground state is determined by comparing 
energies of these states. We find that the antiferromagnetic state is the ground state at 
half filling. When the filling factor is reduced to certain values, the superfluid state takes 
over. Depending on the parameters, different superfluid states, such as d-wave (d), extend 
s-wave (xs), and the mixed-symmetry (xs + i x d) states, can appear. Phase separation 
occurs between superfluid and antiferromagnetic phases. Based on these results, the phase 
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diagram at zero temperature is obtained and its experimental implication is discussed. 



II. EFFECTIVE MODEL 

The Hamiltonian of a single 40 K 87 Rb molecule trapped in a two-dimensional optical lattice 
in the presence of an AC electric field is given by 

Ho = h Hop t + H rot i (1) 

2m 

where the first r.-h.-s. term is the molecule kinetic energy for the center-of-mass motion, and 
the second r.-h.-s. term is the potential energy in the optical lattice. The rotational energy of 
the molecule in the presence of the AC electric field is given by H rot = BL 2 — d-EAc{t), where 
L is the angular momentum operator, B is the rotational constant, d is dipole operator, and 
E j 4c(t) is the AC electric field. The eigenstate of the angular momentum is given by \j, m), 
where j and m are quantum numbers of the total angular momentum and its z-component. 
In this work we consider a linearly polarized AC field in the z direction with amplitude Eac, 
^Acif) = Eac cos(u)i)z. The frequency u is set to couple rotational states |0,0) and |1,0), 
(0, 0|d • E^c(t)|l, 0) = AEac cos (cot) /a/3, where d is the dipole moment. Under the rotating 
wave approximation [7| in the basis of 1 , 0) and exp(— icot)\l, 0), the molecule rotational 
energy is given by 

( Q 

H ro t ~ 



V 



O 2B - Huj 

where f2 = dE^c I v3- In the following we concentrate on the resonance limit to — > 2B / h 
with small electric field, where the two rotational states are nearly degenerate. 
The dipole-dipole interaction is given by 

, \ (1, •(!; 3jd, -e,,) 

v dd [Yi - Tj) 1 _ ,3 ' y z > 

where dj is a dipole at r-j, and is the direction of r^ — rj. In the presence of an resonant AC 
electric field, only rotational states |0, 0) and |1, 0) are important. The matrix elements of the 
dipole operator in this subspace are given by (0, 0|d|0, 0) = (1, 0|d|l, 0) = (0, Old^ll, 0) = 0, 
and (0, 0|rf 2 |l, 0) = d/\/2>. For this two-level system, we can define the pseudospin states 
| |) = (|0,0) + |l,0))/ v / 2 and | \) = (|0,0) - (1,0))/^, so that the dipole matrix elements 
are (f |d| \) = 0, (f |d| f) = dz/y3, and (\, |d| \) = —dz/\/3. In this new basis, the 
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interaction between two dipoles can be conveniently written as a dipolar Ising interaction 
between pseudospins, 

, * 7 o 1 ^ COS Oij /o\ 

V dd {Yi - Tj) = Ad — 7^s iz s jz (3) 

6 I r i r j I 

where Si Z is the z component of the pseudospin s i; and % is the angle between and the z 
direction. For this system, the optical lattice is set in the plane perpendicular to the electric 
field, 9ij = tt/2. 

For a dipolar Fermi gas trapped in a two-dimensional optical lattice in the presence of a 
resonant AC electric field, the effective Hamiltonian is therefore given by 

i 

H = ~ 12 K4a C 3<T + c ja c i<r) + £ 12 J ij s iz s 3* + U 12 n ^ n ih ( 4 ) 

where q ct is the fermion annihilation operator for site % and spin a, rii a = c\ a Ci a , and the 
spin coupling constant between fermions at positions rj and is given by = 4d 2 /(3\rj — 
rj| 3 ). For simplicity we have assumed that in the optical lattice the fermion can only hop 
between nearest neighbors with hopping amplitude t. The onsite-interaction U in Eq. (jlj) 
includes both the s-wave scattering between different spin components and the onsite dipole 
interaction, and it can be tuned either by using Feshbach resonance or by changing the 



vertical confinement 



, [9] . Here we concentrate on the effect of the dipolar Ising interaction 
and set U = in the rest of this paper. 

Since the dipolar Ising interaction is a short-range interaction in two dimension and the 
spin coupling constant is the largest between nearest-neighbor spins, the system tends to 
lower its energy by placing fermions with opposite spins at nearest-neighbor sites, which 
can lead to two different types of ground states. The first one is the antiferromagnetic 
state in which spins align oppositely at nearest neighbors. The other type is the superfluid 
state with exotic symmetries where a fermion attracts another fermion with the opposite 
spin at nearest neighbors to form a Cooper pair. Under the current experimental condition 
[1] , the strength of the dipolar Ising interaction between nearest neighbors is approximately 
J = Ad 2 /(3a 3 ) = 300hHz for the singlet rovibrational state, the hopping amplitude t is 



close to lOOh Hz [ll|] for a typical optical lattice with the strength Vq/Er = 10, where 
Vl(x) = Vosm(k x x), E R = h 2 k 2 ,/2m, and the lattice constant is a = 1/im. In contrast, 
the super-exchange interaction of an ordinary Fermi gas is usually much smaller than the 
hopping amplitude t, which makes it more difficult to produce ordered states. The dipolar 
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Fermi gas has another advantage that the spin fluctuation in Ising systems is much less than 
that in Heisenberg systems, which is helpful to stabilize ordered states. For these reasons, 
the superfluid and antiferromagnetic states may be easier to be created in the dipolar Fermi 
gas than in an ordinary Fermi gas. In the following we focus on the properties of these 
ordered states at zero temperature. 



III. ANTIFERROMAGNETIC STATE 



To study the antiferromagnetic state, we divide the square lattice into two sublattices A 
and B so that nearest neighbors of one sublattice all belong to the other sublattice. The 
antiferromagnetic order is described by the finite average of the pseudospin, (s^) = m for 
sublattice A and (sj Z ) = —m for sublattice B. In the k-space, the mean-field Hamiltonian 
for the antiferromagnetic state is given by 

h af = -^Y.( a U a ^ + b U b U- a U a H- b U b w) 
z k 

+ £ k(4A CT + & L fl ka) + km N ) (5) 

k,cr - 

where N is the total number of lattice sites, the k-summation is over the first Brillouin 
zone of a sublattice, and b^ a are fermion annihilation operators for sublattice A and B, 
£k = — 2t(cos k x a + cos k y a), and a is the lattice constant. The constant h can be regarded 
as the strength of the effective antiferromagnetic interaction, h = — X^o( — l) 4 «^0t — 2. 646 J. 

The mean-field Hamiltonian in Eq. dSJ) can be diagonalized by a standard canonical 
transformation, 

Haf = E( £ k«L«k CT + e£PLM + 2 , (6) 

k,cr - 

where the quasi-particle operators are given by = u^a^ + v^b^, = u^b^i + v^ay^, 
A<t = M k^kt — ^k^kt; A4 = u k a ki ~ v k b ki: an d the transformation coefficients are given by 
u\ = 1 — i>k = [1 + /im/(2ejjr)]/2. The quasi-particles form two bands with energies given by 



^ = ±fk + (^) 2 - (7) 
The magnetic moment m can be determined self-consistently from the equation 

m = 12 ( a kt fl kt - a li a *l) > ( 8 ) 
JV k 



which leads to the equation 

where /i is the chemical potential. Eq. (Q is equivalent to the energy-extreme condition 
BE /dm = 0, and can be solved together with the equation for the fermion number per site 

™ = |B^-^) + ^-4)]- (io) 

The total energy per site is given by 

e af = EMC" - ^) + - 4)} + —■ (ii) 

At half filling n — 1, the r.-h.-s. of Eq. ([9]) is infraredly divergent at m — 0, indicating 
that the antiferromagnetic order always exists for any finite dipolar Ising interaction. The 
antiferromagnetic solution can also be obtained at other filling factors as long as the following 
condition is satisfied, 

^<^E^(A*+|ekl)-f(A*-|ek|)]. (12) 

However, as shown in Fig. Q], the compressibility of the antiferromagnetic state turns nega- 
tive, indicating that the antiferromagnetic state is dynamically unstable beyond half filling. 



IV. SUPERFLUID STATES 

At the superfluid state, two fermions with opposite spins form a pair. In the Hamiltonian 
in Eq. 01]), this pairing interaction is given by — YU^j Jij c lf c it c ]i c jl/ ^- l n the k-space, the 
mean-field Hamiltonian for the superfluid state is given by 

h sf = E( £ k - /^)cL c ka + E( A kC-k+c kt + A k 4 t cL H - A k # k ), (13) 

ktr k 

where the pairing gap is given by A k = £ k , fwgu, g k = (c- k ;c kt ), and 

/kk' = ^ Ji i ex P^ k - k ') • ( r * - r i)]- 

The mean-field Hamiltonian Eq. (fTBl can be diagonalized by Bogoliubov transformation 
as in the BCS theory, 

H SF = N{E SF - /m) + E E k dl a d ka , (14) 

ko- 
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FIG. 1: The energy of the antiferromagnetic state as a function of the filling factor n at J/t = 3. 
The downward curvature shows that the compressibility is negative and antiferromagnetic state is 
unstable beyond half filling. 



where E k = J (e k — /x) 2 + | A k | 2 is the quasiparticle energy, d kCT is the quasi-particle annihi- 
lation operator, and E$f is the ground-state energy per site, 

e sf = ^ - A* - + + tm. (15) 

k k 

The gap A k can be determined from the gap equation 
The fermion number per site is given by 

^SV- (17) 



n 



We numerically solve the gap equation (|T6|) together with Eq. ffTTj) . and find that the gap 
A k can display different symmetries depending on the parameters. The d-wave superfluid 
state appears in the region with 0.5 < n < 1 and J < 51 As shown in Fig. |2fa), the 
d-wave symmetry is characterized by the sign change in the gap at k x ± k y = 0. The 
extended-s-wave (xs) superfluid state may appear for n < 0.5 with the sign change of the 
gap at k x ± k y = ir /a as shown in[2]^b). Between the d-wave and extended s-wave phases, a 
superfluid phase appears with a mixed xs + i x d symmetry in the order-parameter. In this 
phase, a global unitary transformation can be always applied so that the real part of the 
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FIG. 2: The gap A^ versus the x component of the wavevector k x . The dashed and solid lines are 
for Ak at k x = k y and k x + k y = tt/cl respectively, (a) At J/t = 3 and n = 0.656, the gap A^ shows 
the d-wave symmetry, (b) At J/t = 3 and n = 0.29, the gap shows the extended-s-wave symmetry. 
At J/t = 3 and n = 0.483, (c) the real part of the gap shows the extended-s-wave symmetry, and 
(d) the imaginary part of the gap shows d-wave symmetry. 



gap has the extended-s-wave symmetry, and the imaginary part of the gap has the d-wave 
symmetry, as shown in Fig. [2]^c) and|2]^d). Unlike the d-wave or extend-s-wave state, this 
mixed-symmetry state breaks the time-reversal symmetry, and supports the spontaneous 
current 10 ]. 



V. ZERO- TEMPERATURE PHASE DIAGRAM 



By comparing the energies of the antiferromagnetic and superfluid states, we obtain the 
zero-temperature phase diagram of this system, as shown in Fig. [3j At half filling, the 
antiferromagnetic state is the ground state, but it is unstable beyond half filling. When 
J < 5t, the d-wave superfluid state appears at a smaller filling factor. Phase separation 
between the antiferromagnetic and d-wave superfluid states takes place in between. As the 
filling factor further decreases, the d-wave superfluid state evolves into the mixed-symmetry 
(xs + i x d) state. As the filling factor continues to decrease, this mixed-symmetry eventually 
evolves into the extended-s-wave superfluid state. 

When J > 5t, the above picture slightly changes. The d-wave superfluid state no longer 
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FIG. 3: Phase diagram at zero temperature. The antiferromagnetic phase is along the half-filling 
line. The dashed-dotted line marks the boundary between the phase-separation region and the 
superfluid phase. The dashed (solid) line is the phase boundary between the mixed-symmetry and 
the d-wave (extend-s-wave) superfluid phases. 



shows up. As the filling factor decreases from half filling, a superfluid state appears with 
the mixed-symmetry. Phase separation between the antiferromagnetic and mixed-symmetry 
superfluid states occurs in the middle. As the filling factor further decreases, this mixed- 
symmetry evolves into the extended-s-wave superfluid state. 

In the current experimental situation [1], the strength of the dipolar Ising interaction is 
approximately J = 300hHz for the singlet rovibrational state, and the hopping amplitude 
is approximately lOOhHz for a typical optical lattice with lattice constant a = Ifim. 
When the filling factor changes from zero to half filling, extended-s-wave, mixed-symmetry, 
and d-wave superfluid phases, the phase-separation region, and the antiferromagnetic state 
at half filling can all be reached as shown in Fig. |3j The hopping amplitude t can be tuned 
by increasing or decreasing the strength of the optical-lattice potential, and other regions in 
the phase diagram can also be reached. 



VI. DISCUSSION AND CONCLUSION 



So far we have considered only the case for vanishing onsite interaction U — 0. With finite 
U, the phase diagram can be more complicated. If U > and U < J, the antiferromagnetic 
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state is more energetically favorable, and we expect a wider region of the phase separation 
between the antiferromagnetic and superfluid states in the phase diagram. If U ^> J, for 
low energy states there is effectively a constraint of no more than one fermion per site as 
discussed in the t — J — V — W model [3] , which can strongly affect the superfluid states. If 
U < 0, a pure s-wave superfluid state may appear depending on the parameters. All these 
cases are worth to explore in future works. 

Experimentally the onsite interaction leads to losses of the KRb molecules through chem- 
ical reaction, which pose difficulties for creating the antiferromagnetic and superfluid phases 
proposed in this paper. Progress has been made to suppress the reaction rate in a quasi-2D 
;her polar molecules systems, such as KCs, NaK and NaRb, there are 
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no reactive channels 14j , and the loss problem can be avoided. With more experimental 
progresses, it is hopeful for observing the phase diagram shown in Fig. |3j 

In summary, we have shown that in the presence of a resonant AC electric field, the 
dipole-dipole interaction in a dipolar Fermi gas can be turned to a dipolar Ising interaction 
between pseudospins. When trapped in a square optical lattice, the dominant antiferromag- 
netic interaction between nearest neighbors drives the system into different ordered states 
depending on parameters. The antiferromagnetic state is the ground state at half filling. The 
superfluid state appears as the ground state at a smaller filling factor. In the middle is the 
phase-separated region. Depending on the parameters, the superfluid order parameter dis- 
plays different symmetries, such as the d-wave, extended-s-wave, and the mixed symmetry. 
These states mentioned above may be detectable under experimental conditions. 
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